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LIMITING  DETECTION  PERFORMANCE  FOR  RANDOM  SIGNALS 
OF  UNKNOWN  LOCATION,  STRUCTURE,  EXTENT,  AND  STRENGTH 

INTRODUCTION 

This  technical  report  is  the  fourth  in  a  series  of  NUWC 
technical  reports  by  this  author,  covering  the  following  topics: 

(a)  modified  generalized  likelihood  ratio  processors, 

(b)  generalized  likelihood  ratio  processors, 

(c)  power-law  processors,  and 

(d)  optimum  processing, 

respectively.  Topic  (a)  was  analyzed  in  [1],  resulting  in  a 
substantial  compilation  of  receiver  operating  characteristics  for 
the  breakpoint  modification  considered  there.  Topic  (b)  was 
addressed  in  [2],  again  resulting  in  numerous  receiver  operating 
characteristics  that  quantify  the  performance  of  the  modification 
called  the  sum-of-M-largest  processor.  Topic  (c)  was  studied  in 
[3],  resulting  in  a  recommendation  for  a  device  with  power-law 
2.5,  which  performed  remarkably  well,  regardless  of  the  (unknown) 
number  of  bins  occupied  by  signal  and  the  (unknown  but  common) 
signal  level. 

The  overall  goal  of  the  extended  investigation  is  to 
determine  classes  of  processors  that  perform  at  or  near  the 
optimum  level  of  performance,  and  which  can  be  easily  realized 
and  analyzed,  even  in  these  situations  of  scant  knowledge  about 
the  detailed  signal  characteristics.  The  reader  should  be 
familiar  with  the  earlier  material  and  methods  before  undertaking 
the  current  analyses  and  results. 
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PROBLEM  DEFINITION 

The  search  space  consists  of  N  (frequency)  bins,  each 
containing  independent  identically  distributed  noises  of  unit 
power  under  hypothesis  Hq,  signal  absent.  This  situation  is 
presumed  to  be  accomplished  by  an  earlier  normalization 
procedure.  The  number  N  is  under  our  control  and  is  always  a 
known  quantity.  When  signal  is  absent,  the  probability  density 
function  of  each  of  the  bin  output  noises  is  completely  known. 

When  signal  is  present,  hypothesis  H^,  the  quantity  M  is  the 
actual  number  of  bins  occupied  by  the  signal;  this  is  frequently 
an  unknown  parameter.  The  quantity  L  is  the  actual  set  of  bins 
occupied  by  signal,  when  signal  is  present;  for  example,  if 
M  =  4,  then  we  might  have  for  the  occupied  set,  L  =  {2,3,7,29}, 
meaning  that  bins  2,3,7,29  have  signal  in  them.  This  quantity  L 
is  always  unknown  in  our  investigation.  Finally,  the  quantities 
{ Sjn }  are  the  actual  average  signal  powers  in  the  m-th  bin  in 
occupied  set  L,  when  signal  is  present;  these  average  signal 
powers  are  unknown.  We  shall  presume  here  that  all  the  actual 
signal  powers  per  bin  are  equal  to  a  common  (unknown)  value  S  in 
the  occupied  set  of  bins,  L,  and  zero  elsewhere. 
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Observe  that  the  actual  signal  power  per  bin,  S,  can  also  be 
interpreted  as  the  actual  signal-to-noise  power  ratio  per  bin, 
since  the  noise  power  per  bin  has  been  normalized  at  unity. 


JOINT  PROBABILITY  DENSITY  FUNCTIONS  OF  THE  SET  OF  OBSERVATIONS 

The  joint  probability  density  function  governing  the 
complete  observation  {xnJ  under  hypothesis  Hq  follows  from 
(1)  and  the  statistical  independence  of  {xn)  as 

N 

Pq^,  . . .  ,uN)  =  fT{exp(-un)  J  .  (4; 

n=l 

Under  hypothesis  H1,  the  signal  can  land  in  any  set  of  M 
disjoint  bins,  out  of  the  total  of  N  search  bins.  This  results 
in  a  total  number  of  possibilities  K  =  (N|M),  where  the  quantity 
in  parentheses  is  the  binomial  coefficient.  Without  any  apriori 
information  about  signal  location,  it  is  presumed  that  each  set 
occurs  with  equal  probability  1/K.  Thus,  there  are  K  possible 
occupancy  sets,  each  of  size  M,  namely  sets  {L^}  for  1  <  k  1  K. 
The  joint  probability  density  function  governing  the  complete 
observation  {xnJ  is  therefore 

K  fl  , 

■  EH  K  TTla  exp(-aun)  }  TTlexp(-un)}  ,  (5) 

k=l  neLk  n^Lk 

where  we  used  ( 1 )  and  ( 2 ) . 
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SIGNALS  WITH  PARTIAL  LOCATION  INFORMATION 

In  this  section/  the  total  number  of  bins,  M,  occupied  by  a 
signal  (when  present)  and  the  (equal)  average  signal  powers  per 
bin,  S,  are  presumed  known  to  the  optimum  processor.  However, 
the  locations  of  the  M  occupied  bins,  in  the  total  of  N  search 
bins,  are  only  partially  known  to  the  optimum  processor.  Instead 
of  the  M  occupied  bins  being  allowed  to  occur  anywhere  in  the 
range  of  N  bins  (which  would  mean  a  number  of  possibilities  equal 
to  binomial  coefficient  K  =  (N|M)),  we  presume  in  this  section 
that  a  lesser  number  of  possibilities,  Kp,  for  the  occupied  bin 
patterns  can  occur.  Further,  each  of  these  Kp  allowed  patterns 
occurs  with  equal  probability  1/Kp.  In  particular,  each  location 
pattern  Lk,  1  <  k  £  Kp,  is  a  set  of  size  M,  with  a  structure 
yet  to  be  specified.  For  example,  if  M  *  4,  we  might  have 
L7  =  {5,6,13,19},  meaning  that,  out  of  the  total  of  N  search 
bins,  the  particular  bins  5,6,13,19  are  occupied  by  signal,  for 
pattern  Lj . 
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IIZ  exp(n  xk)  <  v  •  (10) 

k- 1 

The  likelihood  ratio  test  indicates  that  we  must  compute  all  the 
K  possible  allowed  linear  sums  {X,  }  in  (9),  weight  each  sum  by 

p  K 

common  value  w,  exponentiate,  and  sum  over  all  the  possible  K 

XT 

sets.  However,  Kp  here  may  be  smaller  than  the  earlier  value 
K  =  (N|M)  encountered  in  (5),  due  to  fewer  allowed  set  patterns 
here. 

There  has  been  no  restriction  on  the  form  of  the  Kp  sets  {L^l 
that  describe  the  allowed  occupancy  patterns  for  the  M  occupied 
bins,  except  that  each  set  is  equally  likely  and  is  of  size  M. 
If  these  Kp  sets  are  carefully  chosen,  the  difficult  and  lengthy 
calculation  prescribed  in  (9)  -  (10)  for  the  optimum  processor 
can  be  greatly  simplified,  leading  to  a  workable  procedure 
capable  of  simulation  in  a  reasonable  amount  of  time,  for  some 
combinations  of  parameter  values. 

If  we  define  random  variables 


en  =  exp(w  xn)  for  1  £  n  £  N  , 


(11) 


the  likelihood  ratio  test  in  (10)  can  be  expressed  as 


if. 

k=l  ncL, 


'n  < 


v 


(12) 


Each  of  the  Kp  products  is  of  size  M.  This  form  will  prove 
useful  in  later  manipulations,  where  we  will  manage  to  convert 
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GENERALIZED  LIKELIHOOD  RATIO  PROCESSOR 

In  this  subsection,  we  presume  that  M,  the  number  of  occupied 
signal  bins  (when  present)  is  known,  but  that  the  actual  average 
signal  strength  S  and  the  particular  occupied  bin  locations  L  are 
not  known.  Therefore,  we  hypothesize  average  signal  strength 
S  k  0  and  occupied  set  Lk,  where  we  let  1  i  k  <  K^,  as  in  the 
previous  subsection.  The  size  of  every  set  Lk  is  M.  Letting 
a  =  1/(1+S)  <  1,  the  governing  joint  probability  density 
function  under  hypothesis  is,  for  the  hypothesized  S  and  L^, 

P1(u1,...,uN)  =  TTla  exp(-aun)  }  TT l exP ( "un )  1  •  (13) 

neLj^  n/L^ 

For  observation  fx  },  define  the  random  variables 

n 

N 

Z  =  C  xn  ,  Xk5C*n  forl<k<K  .  (14) 

n=l  neLk  * 

Then,  the  joint  density  in  (13)  takes  on  value 

P1  s  Pi(xi^---/xN)  *  a“  exp(-aXR)  exp[-(Z-Xk)]  = 

=  exp(-I)  a-  exp[(l-a)Xk)  ,  (15) 

which  depends  on  a  and  k.  Holding  hypothesized  value  k  fixed, 
the  choice  of  continuous  parameter  a  that  maximizes  in  (15)  is 
random  variable 

a  f»/Xk  ifXk>M, 
l  1  if  xk  <  m] 
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Using  the  monotonicity  of  functions  exp  and  g  in  (19),  the 
generalized  likelihood  ratio  test  is  therefore  simply 

max  =  max |x^ , . . . ,XR  j  *  v  .  (22) 

Physically,  this  test  for  signal  presence  makes  a  great  deal 
of  sense.  First,  each  variable  X^  in  (14)  is  obtained  by 
performing  as  much  linear  summation  over  the  observation  {xnJ  as 
possible,  namely  over  the  particular  M  random  variables  that 
constitute  each  set  Lk,  thereby  taking  advantage  of  the  known 
limited  signal  structure  information.  Then,  the  largest  sum,  out 
of  the  total  of  possibilities,  is  compared  with  a  threshold 
for  a  decision  about  signal  presence  or  absence.  In  this 
generalized  likelihood  ratio  test  (22),  there  is  no  need  to  know 
the  actual  average  signal  power  S  or  the  actual  locations  L  of 
the  occupied  bins.  However,  M  must  be  known  and,  along  with  the 
partial  signal  location  information  {L^},  is  used  to  form  {X^J  in 
(14). 

This  test  in  (22)  is  not  identical  to  the  optimum  likelihood 
ratio  test  (10)  or  (12).  However,  with  careful  selection  of  Kp 
and  the  attendant  allowed  sets  { ,  the  generalized  likelihood 
ratio  test  in  (22)  is  often  easy  to  realize  in  practice,  and  can 
be  simulated  in  a  reasonable  amount  of  time  in  order  to 
accurately  ascertain  its  performance,  whereas  optimum  test  (12) 
frequently  cannot.  Also,  the  optimum  test  (12)  requires 
knowledge  of  the  actual  average  signal  power  per  bin,  S. 
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SIGNALS  WITH  ZERO  LOCATION  INFORMATION 

Here,  the  number  of  bins,  M,  occupied  by  signal  and  the 
(equal)  average  signal  powers  per  bin,  S,  are  presumed  known  to 
the  optimum  processor,  but  the  locations  of  the  occupied  bins  are 
completely  unknown.  The  signal  can  land  in  any  set  of  M  disjoint 
bins,  each  set  occurring  with  equal  probability  1/K,  where  K  is 
the  binomial  coefficient  (N|M).  There  are  now  K  occupancy  sets 
{L,},  each  of  size  M. 


OPTIMUM  PROCESSOR 


Under  hypothesis  HQ,  the  governing  probability  density 


function  is 


P0(u1,...,uN)  =  Tn_{exp(-un)  J 

n— 1 


(24) 


Under  hypothesis  ,  letting  a  «  1/(1  +  53) ,  we  have,  instead, 
density 


K  fl  , 

p1(u1, . . .  ,uN)  =  Til  [k  TTla  exp(-aun)}  TT{exp(-un)  }  j  .  (25) 

k— 1  neL,  n/L, 


The  likelihood  ratio  for  observation  {x„l  is 

n 


LR  = 


Pl(Xi,---,* N)  1  M  K 

d~7x - FT  =  k  —  ■  1  T { exP ( x  [  1  -  a])} 

P0(x1,...,xN)  k  k=1  neL^  n 
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GENERALIZED  LIKELIHOOD  RATIO  PROCESSOR 

In  the  case  where  M  is  known,  but  average  signal  power  S  is 
not,  we  can  appeal  directly  to  results  (22)  and  (14)  for  the 
generalized  likelihood  ratio  test  and  interpret  the  parameters 
differently.  Namely,  we  obtain  the  test 

max[x1, . . . ,XK)  <  v  ,  (31) 

where 

Xk  =  E  xn  for  1  <  k  S  K  =  [JJ]  .  (32) 

neLk  v-' 

Notice  that  K  is  now  the  very  large  integer,  (N|M),  required  when 
there  is  no  information  on  the  signal  location. 

On  the  other  hand,  if  both  M  and  S  are  unknown,  we  resort  to 
(23)  and  (14),  which  gives  us  the  generalized  likelihood  ratio 
test 

p  ' 

max  M  gfmax  X,/m]  ^  v  ,  (33) 

'  k  K  ;J  ^ 

but  now  where 

Xk  =  E  xn  for  1  i  k  i  K(M)  =  (JJ)  .  (34) 

neLk  '  v  ' 

Here,  K(M)  =  (N|M)  is  a  function  of  the  hypothesized  number,  M, 
of  occupied  bins . 
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SIGNALS  WITH  SLIDING  LOCATION  INFORMATION 

Here,  we  again  presume  M  and  S  are  known,  but  now  we  restrict 
the  occupied  bin  structure  to  be  M  contiguous  bins  of  unknown 
(sliding)  location.  There  are  now  only 

K  =  N  -  M  +  1  (35) 

c  —  '  ' 

possible  locations  for  the  contiguous  band,  which  are  assumed 
to  be  equally  probable,  apriori.  Let  the  k-th  set  now  be 

Lk  =  {k,k+l, . . . ,k+M-l }  for  1  <  k  S  Kc  .  (36) 

Physically,  this  corresponds  to  a  signal  of  known  bandwidth,  but 
of  unknown  center  frequency. 
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likelihood  ratio  test  (where  £5  is  no  longer  needed) 

Kc  ,  , 

<  < v  '  <41> 

with  Xk  given  by  the  sliding  sum  (39)  of  M  linear  terms.  This 
approximate  likelihood  ratio  test  (41)  linearly  accumulates  the 
data  {xnJ  as  much  as  possible,  to  yield  the  partial  sums  {X^}  in 
(39),  but  then  adds  up  quadratic  versions  of  these  latter  sums. 
Implementation  of  test  (41)  is  a  relatively  simple  task,  because 
K  is  generally  a  manageable  number  here.  A  worthwhile 

v  * 

alternative  is 

Kc 

C  xj  >  v  ,  (42) 

but  where  the  best  choice  of  power  p  is  not  obvious. 

It  should  be  noted  that  the  approximate  likelihood  ratio 
tests  (41)  and  (42)  use  a  combination  of  all  the  partial  sums 
{XkJ  in  making  their  decision  about  signal  presence  or  absence; 
however,  the  power-law  operation  accents  the  larger  partial  sums 
at  the  expense  of  the  smaller  ones.  Also,  tests  (41)  and  (42) 
require  very  little  computation  and  can  be  implemented  without 
knowledge  of  average  signal  power  level  S;  of  course,  the  number 
M  of  occupied  bins  must  be  known  in  order  to  evaluate  (39). 
Simulation  comparisons  of  (40)  -  (42)  would  have  to  be  conducted 
in  order  to  ascertain  the  exact  degradations  in  performance 
caused  by  adoption  of  the  suboptimum  processors. 
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OPTIMUM  BANDED  PROCESSOR 

We  will  now  apply  the  general  result  for  the  likelihood  ratio 
test  in  (10)  and  (12)  to  the  particular  example  where  the  allowed 
Kp  patterns  for  the  M  occupied  bins  are  formed  from  a  banded 
signal  structure  in  frequency  space.  In  particular,  we  divide 
the  total  search  space  of  N  bins  into  M  disjoint  bands,  each 
containing  I  =  N/M  bins.  It  is  presumed  that  I  is  an  integer. 

We  now  consider  a  signal,  structured  so  that  one,  and  only 
one  signal,  can  lie  in  each  band  of  I  bins.  This  does  not 
correspond  to  the  practical  physical  situation,  where  the  M 
signals  can  lie  anywhere;  rather,  this  signal  model  is  concocted 
solely  for  purposes  of  simplifying  likelihood  ratio  test  (12). 

It  is  recognized  that  the  performance  of  the  corresponding 
optimum  banded  processor  (OPT-B)  will  be  better  than  the  optimum 
processor  with  zero  location  information  (OPT-Z),  namely  (28). 
Therefore,  the  performance  of  the  OPT-B  will  furnish  a  bound  on 
performance  for  any  practical  processor;  hopefully,  it  will  also 
turn  out  to  be  a  useful  and  tight  bound. 

By  restricting  the  M  signal  components  to  lie  in  disjoint 
equal-size  bands,  and  then  giving  this  information  to  the  OPT-B, 
some  partial  location  information  has  been  furnished.  This 
information  is  in  addition  to  the  values  of  M  and  S,  which  are 
also  presumed  known  to  the  OPT-B.  However,  to  compensate  as  much 
as  possible  for  this  additional  information,  we  require  each 
individual  signal  component  to  be  equally  likely  to  occupy  any  of 
the  I  =  N/M  bins  in  each  band.  This  maximizes  the  confusion  for 
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optimally.  Since  optimum  is  defined  in  terms  of  maximum  P^  for  a 
specified  P^,  better  performance  is  reflected  by  a  curve  being 
nearer  the  top  of  the  figure.  On  the  other  hand,  the  bottommost 
curve  for  the  best  power-law  device  (best  v)  must  always  lie 
below  the  OPT-Z,  because  the  power-law  device  has  no  knowledge  of 
M  or  S,  nor  is  it  optimum  in  any  sense. 

We  will  demonstrate,  by  means  of  several  numerical  examples, 
that  the  top  and  bottom  curves  in  figure  1  lie  virtually  on  top 
of  each  other.  That  is,  the  best  power-law  device  performs  so 
well,  and  the  upper  bound  furnished  by  the  OPT-B  is  so  tight, 
that  the  region  between  these  two  is  a  very  narrow  ribbon.  In 
fact,  it  will  be  shown  that  if  the  signal  power  S,  for  the  best 
power-law  device,  is  increased  by  0.1  dB,  then  these  two 
corresponding  receiver  operating  characteristics  overlap.  Thus, 
we  will  have  succeeded  in  pinching  the  unknown  receiver  operating 
characteristic  of  the  OPT-Z  within  a  very  tight  ribbon,  and  will 
have  effectively  determined  its  performance. 

We  will  also  have  accomplished  two  other  worthwhile  goals. 
First,  we  will  have  found  a  practical  device,  namely  the 
power-law  processor,  which  performs  virtually  at  the  optimum 
level  of  performance  in  this  environment.  Second,  we  will  have 
established  such  a  tight  upper  bound  on  performance  (by  means  of 
the  OPT-B)  that  there  is  no  further  need  to  search  for  tighter 
bounds  or  for  the  exact  performance  of  the  OPT-Z.  A  discrepancy 
of  0.1  dB  is  so  small  that  further  significant  effort  on  tighter 
bounds  is  hard  to  justify. 
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where  en  =  exp(w  xn),  w  =  S/(l  +  S) .  Notice  that  the  sum  of 
Kp  =  16  products  in  (47)  has  been  condensed  into  a  single  product 
of  sums  in  (48).  This  collapsing  feature  occurs  as  a  direct 
result  of  the  banded  signal  structure.  Inspection  of  the  final 
form  in  (48)  reveals  that  it  allows  for  interaction  of  every 
component  en/  1  <  n  £  4,  in  the  first  band,  with  every  component 
eR,  5  £  n  £  8,  in  the  second  band,  but  no  interaction  within 
either  band.  This  is  the  mathematical  representation  of  the 
physical  separation  of  signals  into  disjoint  bands. 

This  collapsing  behavior  of  (12)  holds  in  general  for  the 
banded  signal  structure.  Namely,  the  likelihood  ratio  test  in 
(12)  can  be  generally  expressed  exactly  in  this  special  case  as 

p  p  •  •  •  p  4>  t  •  •  4>  p  p  •  •  •  p  r: 

el  eI+l  eN— I+l  T  T  eI  e2I  eN 

=  (ej  +  •••  +  ej|(eI+1  +  •••  +  ®2I  >  ’ ' '  "tj-i+I  +  +  V  "  <49> 

■  n  ■  n  j^r exp  )  < v '  1 50  ’ 

where  there  are  I  terms  in  each  of  the  M  sums.  Observe  that 
expansion  of  (49)  would  yield  I—  =  Kp  terms,  which  is  just  the 
number  of  terms  in  form  (12).  However,  compact  form  (49) 
requires  only  the  evaluation  of  N  exponentials  {enJ  and  M-l 
products,  whereas  form  (12)  requires  N  exponentials  and  Kp  (M-l) 
products,  a  number  which  is  intolerably  large  when  N  is  large. 
Result  (10)  for  the  likelihood  ratio  test  is  not  a  useful 
alternative  because  it  requires  K  exponentials  and  summations. 
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Furthermore,  whereas  we  are  unaware  of  any  method  to  cut 
down  on  the  number  of  calculations,  K,  required  for  the  OPT-Z 
likelihood  ratio  test  (28),  the  OPT-B  likelihood  ratio  test  (12) 
can  be  reduced  to  form  (49),  which  is  very  easy  to  simulate  in  a 
reasonable  amount  of  computer  time. 

As  another  example,  for  N  =  1024,  M  =  4,  there  follows 
K  =  4.6E10  and  =  4.3E9.  Now,  the  number  of  alternatives  has 
been  reduced  by  a  factor  of  10,  and  the  adequacy  of  the  tightness 
of  the  performance  bound  for  the  OPT-B  becomes  even  more 
questionable.  On  the  other  hand,  there  are  still  =  4.3E9 
alternatives  left  for  the  banded  signal  structure,  which 
indicates  that  a  great  deal  of  uncertainty  still  remains  for 
the  OPT-B  to  consider. 

Finally,  for  N  =  1024  and  M  =  512,  we  find  K  =  4.5E306 
and  K  =  1.3E154.  Here,  the  number  of  alternatives  K  is 

p  p 

approximately  the  square  root  of  K;  that  is,  Kp  is  reduced  by  a 
152 

factor  of  10  relative  to  K.  On  the  other  hand,  there  still 
remains  a  tremendous  number  of  alternatives  for  the  OPT-B  to 
consider,  namely  on  the  order  of  10*^,  which  may  be  sufficient 
to  confuse  the  OPT-B  enough  to  yield  a  tight  performance  bound. 
This  question  can  only  be  resolved  through  numerical 
investigation  of  the  receiver  operating  characteristics. 

(For  N  =  M  =  1024,  we  have  K  =  Kp  =  1,  and  the  OPT-B  and  OPT-Z 
processors  become  identical,  namely  the  energy  detector.) 

All  the  above  results  hold  for  the  present  case  of  equal 
signal  powers  S  per  bin.  However,  when  the  signal  powers  per  bin 
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RECEIVER  OPERATING  CHARACTERISTICS  FOR  OPTIMUM  BANDED  PROCESSOR 

We  now  determine  the  performance  of  the  optimum  banded 
processor,  as  characterized  by  (49)  or  (50)  in  general,  for  any 
N  and  M  (with  N/M  integer).  The  weight  w  in  likelihood  ratio 
test  (50)  is  given  by  w  =  S/(l  +  S) ,  which  depends  on  the 
particular  signal  power  case  under  investigation.  This  means 
that,  in  simulating  (50),  the  false  alarm  probability  Pf  (as  well 
as  the  detection  probability  P^)  must  be  rerun  for  each  different 
value  of  signal  power  S  of  interest.  (This  is  distinct  from  the 
power-law  processor,  for  example,  where  one  curve  of  P^  could  be 
plotted  versus  numerous  P^  curves;  also,  P^  was  available 
analytically  for  some  power-law  devices,  which  circumvented  the 
need  for  simulating  Pf  at  all  in  those  cases.) 

This  additional  simulation  expense  for  the  optimum  banded 
processor  leads  us  to  adopt  the  following  procedure.  We  first 
refer  to  the  performance  of  the  best  power-law  device  [3]  to 
determine  what  signal  power  S  was  necessary  there  in  order  to 

hr 

operate  in  the  neighborhood  of  the  low-quality  operating  point 
-3 

Pj  =  10  and  P^  =  0.5.  Then,  we  choose  a  pair  of  signal  powers 
that  bracket  S  :  that  is,  we  take  S„  <  £  S.  and  we  simulate 

(50)  for  both  Sa  as  well  as  S^.  This  gives  us  a  pair  of  receiver 
operating  characteristics  for  the  optimum  banded  processor, 
namely  two  curves  of  P^  versus  Pf,  from  which  we  can  interpolate 
in  order  to  estimate  the  signal  power  needed  to  operate  at  the 
low-quality  operating  point.  This  entire  procedure  was  repeated 
for  the  high-quality  operating  point  P^  =  10-^  and  P^  =  0.9. 


31 


TR  10839 


A  table  of  the  signal  powers  S  in  decibels  required  for  the 
optimum  banded  processor  has  been  extracted  from  these  figures; 
it  is  presented  under  the  column  labeled  OPT-B  in  table  1  below. 
Also  listed  are  the  corresponding  signal  powers  required  for  the 
best  power-law  processor,  under  the  heading  BEST  PLP;  the  last 
column  gives  the  actual  best  value  of  power  v  to  use.  These 
latter  values  come  from  [3;  page  28,  table  1].  The  largest 
discrepancy  in  table  1  between  OPT-B  and  BEST  PLP  is  0.17  dB; 
however,  the  coarseness  of  the  search  that  we  conducted  in  v 
values  (namely  1,  2,  2.5,  3,  ®)  leads  us  to  believe  that  the 
largest  discrepancy  would  be  closer  to  0.1  dB  if  intermediate 
values  of  \>  were  investigated.  Thus,  we  claim  that  the  best 
power-law  processor  is  within  approximately  0.1  dB  of  the 
absolute  limit  of  performance  in  this  environment. 

Table  1.  Required  S(dB)  for  N  =  1024,  P^  =  10-3,  P^  =  0.5 


M 

OPT-B 

BEST  PLP 

V 

1 

12.75 

12.77 

CO 

2 

10.03 

10.11 

00 

4 

7.88 

8.05 

3 

8 

6.04 

6.2 

3 

16 

4.35 

4.4 

3 

32 

2.62 

2.7 

2.5 

64 

0.71 

0.75 

2.5 

128 

-1.54 

-1.4 

2 

256 

-4.16 

-4.0 

2 

512 

-7.03 

-7.00 

1 

1024 

-10.01 

-10.01 

1 
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The  largest  discrepancy  in  table  2  between  OPT-B  and  BEST  PLP 
is  0.3  dB;  however,  the  coarseness  of  the  search  in  v  and  the 
difficulty  of  simulating  extremely  large  trial  runs  for  small 
probabilities  in  the  10-^  range  makes  it  difficult  to  accurately 
estimate  the  difference  in  performance  between  the  optimum  banded 
processor  and  the  best  power-law  processor.  A  reasonable 
estimate  appears  to  be  about  0.2  dB. 

The  difficulty  of  obtaining  accurate  receiver  operating 
characteristics  for  very  small  false  alarm  probability  values  is 
due  to  the  extremely  large  number  of  trials  required.  As  an 
alternative,  an  analytical  approach  for  determining  the  false 
alarm  probability  for  the  optimum  banded  processor  is  outlined  in 
appendix  A.  However,  it  requires  two  very  accurate  numerical 
integrations,  coupled  with  two  fast  Fourier  transforms  of  rather 
large  size.  This  procedure  was  not  employed  here. 

In  figure  22,  the  total  signal  power,  M  S  in  decibels, 
required  to  achieve  the  low-quality  operating  point,  is  plotted 
versus  M  for  the  power-law  processors  v  =  1,  2,  2.5,  3,  when 
N  =  1024.  This  result  comes  from  [3;  page  29],  The  OPT-B 
absolute  bound  in  table  1  above  is  superposed  as  black  dots  on 
this  figure.  As  already  anticipated,  these  results  indicate  that 
there  is  always  some  power-law  processor  that  performs  very  close 
to  optimum. 

A  similar  result  for  the  high-quality  operating  point  is 
given  in  figure  23;  the  source  of  these  results  is  [3;  page  30] 
and  the  OPT-B  absolute  bound  in  table  2  above.  The  danger  of 
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ErQ  ^5  E-4  f-3  ?f  .1  .5- 

Figure  3.  ROC  for  OPT-B,  M  =  2,  N  =  1024,  Pd  ~  0.5 
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Figure  7.  ROC  for  OPT-B,  M 
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Figure  11.  ROC  for  OPT-B,  M  =  512,  N  ■  1024,  P 
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Figure  19.  ROC  for  OPT-B,  M  =  128,  N 
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COMPARISON  WITH  TWO  OTHER  PROCESSORS 

MODIFIED  GENERALIZED  LIKELIHOOD  RATIO  PROCESSOR 

The  modified  generalized  likelihood  ratio  processor  for 
search  size  N  =  1024  was  investigated  in  [1]/  resulting  iri  36 
receiver  operating  characteristics  for  M  =  8,16,32,64,128,256, 
and  breakpoints  xQ  =  1,3,5,7,9,11;  see  [1;  figures  8  -  43], 

The  remaining  36  receiver  operating  characteristics,  for 
M  =  1,2,3,4,512,1024  and  the  same  6  breakpoint  values,  are 
presented  here  in  appendix  B.  A  complete  set  of  characteristics 
for  xQ  =  6  and  M  =  1,2,3,4,8,16,32,64,128,256,512,1024  has  also 
been  added.  All  these  results  were  obtained  by  simulation;  the 
number  of  trials  utilized  is  indicated  on  each  figure. 

From  this  totality  of  results,  the  signal  power  S  required  to 
achieve  specified  false  alarm  and  detection  probabilities  can  be 
extracted.  These  results  are  presented  in  figure  24  for  the 
low-quality  operating  point.  They  indicate  that  the  best  single 
compromise  breakpoint  xQ  for  the  modified  generalized  likelihood 
ratio  processor,  in  order  to  cover  all  M  values  from  1  to  1024, 
is  =  5.6;  the  maximum  losses  at  M  =  1  and  1024  are  3.7  dB 
relative  to  the  OPT-B.  This  is  considerably  poorer  than  the 
power-law  processor  with  power  v  =  2.4,  for  which  the  maximum 
loss  is  only  1.2  dB,  regardless  of  the  value  of  M;  see  figure  22. 


49 


TR  10839 


LINEAR  DEVICE  WITH  DEAD  ZONE 


The  power-law  processor  was  derived  as  an  approximation  to 
the  optimum  processor,  while  the  modified  generalized  likelihood 
ratio  processor  was  the  result  of  an  ad  hoc  change  incorporated 
to  circumvent  weak  signal  estimates.  These  approaches  are 
adopted,  of  necessity,  in  order  to  obtain  realistic  practical 
processors.  However,  they  leave  open  the  possibility  of  yet 
other  suboptimum  approaches  that  may  yield  equally  good  practical 
techniques  for  detection  of  weak  signals  of  unknown  character. 

In  this  subsection,  we  consider  a  processor  that  comes  about 
as  a  modification  to  the  likelihood  ratio  test.  The  derivation 
of  this  processor  (as  well  as  its  performance  in  terms  of  and 
Pd)  is  given  in  appendix  C.  The  end  result  is  the  following 
modified  likelihood  ratio  test,  (C-7),  for  signal  presence: 


N  N 

| 

fx  -  b 

for 

x  >  b 

IZ  h(*n)  <  V  , 

n=i 

h(x)  = 

i  o 

for 

x  <  b 

b  £  0  .  (53) 


The  component  device  h(x)  is  linear  for  x  >  b,  but  has  zero 
response  for  x  <  b.  This  device  completely  suppresses  weak  data 
values  xn,  but  allows  the  stronger  data  values  to  contribute 
linearly.  Selection  of  the  best  choice  of  breakpoint  b  requires 
detailed  computation  of  P^  and  P^  for  numerous  values  of  M  and  S. 

A  closed  form  for  false  alarm  probability  P^  is  given  in 
(C-13).  It  was  used  to  generate  the  series  of  curves  for  P^ 
versus  threshold  v  in  figure  25.  On  the  other  hand,  for  the 
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signal  present  hypothesis  ,  a  closed  form  for  the 
characteristic  function  of  the  decision  variable  is  given  by 
(C-14).  This  latter  result  was  not  used  here  for  determining  P^; 
rather,  P^  was  found  by  simulation,  since  its  most  useful  range 
of  values  is  easily  estimated. 

A  series  of  receiver  operating  characteristics  were 
determined  for  processor  (53),  but  are  not  presented  here,  for 
the  sake  of  brevity.  Tables  of  the  signal  powers  per  bin,  S, 
required  for  this  modified  likelihood  ratio  (MLR)  processor  (53) 
were  extracted  from  these  figures  for  both  the  low-quality  and 
high-quality  operating  points;  they  are  presented  under  the  third 
column  labeled  MLR  in  tables  3  and  4  below.  The  last  column 
gives  the  corresponding  best  value  of  breakpoint  b  to  use.  Also 
listed  are  the  corresponding  signal  powers  required  for  the 
OPT-B.  The  largest  discrepancies  in  tables  3  and  4  between  the 
OPT-B  and  MLR  processors  are  approximately  0.35  dB,  when  the  best 
breakpoint  b  is  used. 

This  modified  likelihood  ratio  processor  needs  to  know  signal 
size  M,  in  order  to  make  the  best  selection  of  breakpoint  b;  this 
unrealistic  condition  is  consistent  with  all  the  processors 
encountered.  However,  we  are  unable  to  make  a  decision  on  the 
best  single  compromise  breakpoint  b  to  use  for  this  MLR 
processor,  because  the  large  number  of  receiver  operating 
characteristics  that  would  be  required  to  cover  all  the  potential 
b  and  M  values  were  not  run.  The  analytic  results  in  appendix  C 
would  be  useful  in  this  regard. 
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SUMMARY 

The  limiting  detection  performance  for  random  signals  of 
unknown  location,  structure,  extent,  and  strength  has  been  very 
accurately  ascertained  through  use  of  a  "super"  optimum  processor 
called  the  optimum  banded  processor.  This  processor  presumes 
knowledge  of  the  number  M  of  occupied  bins,  the  signal  power  per 
bin  S,  and  some  partial  knowledge  about  the  locations  of  the 
occupied  bins.  This  unrealistic  set  of  assumptions  leads  to  an 
upper  bound  on  attainable  detectability  performance  of  any 
practical  processor;  however,  the  bound  turns  out  to  be  so  tight 
that  searching  for  further  improvement  can  yield  only  negligible 
gains,  of  the  order  of  0.1  dB. 

The  tightness  of  the  upper  bound  is  established  by  comparing 
it  with  the  performance  of  the  power-law  processor  using  its  best 
power  v.  It  is  found  that  there  is  always  some  power-law 
processor  (some  value  of  \>)  that  performs  within  0.1  dB  of  the 
absolute  optimum  level.  Strictly,  this  conclusion  is  drawn  only 

_3 

for  the  numerical  example  of  N  =  1024,  P^  =  10  ,  =  0.5.  This 

achievement  for  the  power-law  processor  is  striking,  in  light  of 
the  fact  that  it  knows  and  uses  nothing  about  size  M,  signal 
power  S,  nor  anything  about  the  signal  locations.  However,  when 
it  is  necessary  to  employ  a  single  compromise  value  of  power  v  to 
operate  in  the  entire  range  of  possibilities  from  M  =  1  to 
M  =  N  =  1024,  the  minimax  loss  is  1.2  dB,  occurring  for  power 
selection  v  =  2.4. 


55 


TR  10839 


APPENDIX  A.  FALSE  ALARM  PROBABILITY  FOR  OPTIMUM  BANDED  PROCESSOR 

We  are  interested  in  obtaining  the  false  alarm  probability 
for  the  processor  described  in  (49)  and  (50),  where 

s 

en  -  «xp(w  *n>  -  «  =  fTs  '  “  =  1  +  w  =  2  +  §  '  1  =  5  • 


We  define  the  auxiliary  random  variables 


s  =  >  '  eT  _  .  ,  y  =  ln(s  )  for  1  <  m  <  M  ,  ( A-2 ) 

m  j-'i  rol-I+j  -*m  v  m'  —  '  ' 


z  =  y,  + 


yl  +  ' ' '  +  yM  • 


(A-3) 


The  random  variable  sm  comes  directly  from  (50),  while  z  is  the 
logarithm  of  the  left-hand  side  of  (49).  Since  random  variable 
xn  is  never  negative  here,  we  always  have  en  >  1,  sm  1  I, 
ym  £  ln(  I ) ,  z  >  M  ln(  I ) .  Also,  w  1  2. 

The  fundamental  calculation  is  the  characteristic  function  of 
en,  which  for  £  >  0,  can  be  developed  under  Hq  as 


fe(S)  s  I  i?;enj  =  expliS  exp(wxn)l  = 


=  J  du  exp(-  u  +  i£  exp(wu))  =  (w-1)  J  dt  t  w  exp(iSt)  =  (A-4) 


(w-1 ) ( -i£ ) w  1  J  dy  y-w  exp(-y)  =  (w-1 ) ( -i£ ) W_1  J  dy  y~w  exp(-y) 

-U  -U 


=  (w-1)  (-U)W_1  [T(l-w)  -  Y(  l-w,-i£)  ]  = 
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If  fe(S)  is  evaluated  via  (A-5)  and  (A-6)  at  increment  AJ^, 
then  fs(?)  becomes  available  at  the  same  increment.  Use  of  these 
samples  in  (A-7)  yields  an  approximation  p  (u)  with  period 

which  must  be  large  enough  to  encompass  the  full  extent 
of  ps(u).  If  an  FFT  (fast  Fourier  transform)  of  size  is  used 
to  carry  out  (A-7),  the  increment  for  the  argument  of  p  (u)  is 
Au^  =  2n/(Ni  A£^);  since  Au^  must  be  small,  this  will  require  a 
large  value  for  . 

If  fy(S)  in  (A-8)  is  evaluated  at  increment  A£2,  then  fz(S) 
is  available  at  the  same  increment.  When  these  samples  are  used 
in  the  FFT  required  to  determine  the  density  p  (u)  of  z,  the 

Z 

approximation  Pz(u)  has  period  2n/AE;2,  which  must  be  large.  If 
this  FFT  is  taken  at  size  N2,  the  increment  for  the  argument  of 
pz(u)  is  Au2  =  2n/(N2  AS2) ,  which  must  be  small  enough  to  track 
the  variations. 

The  fundamental  parameters  in  this  investigation  are  N,  M, 
and  S,  where  I  =  N/M  must  be  an  integer.  Two  different  numerical 
integrations  are  needed  to  carry  out  (A-6)  and  (A-8).  Also,  two 
FFTs  are  required:  one  for  p  (u)  from  f  (£)  in  (A-7)  and  the 

W  b 

other  for  p_(u)  from  f „(Z) •  This  intensive  numerical  procedure 
z  z 

will  require  large  values  for  FFT  sizes  and  N2. 
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APPENDIX  B.  RECEIVER  OPERATING  CHARACTERISTICS 
FOR  MODIFIED  GENERALIZED  LIKELIHOOD  RATIO  PROCESSOR 


The  modified  generalized  likelihood  ratio  ( MGLR )  processor 
for  search  size  N  =  1024  was  investigated  in  [1],  resulting  in  36 
receiver  operating  characteristics  for  M  =  8,16,32,64,128,256, 
and  breakpoints  xD  =  1,3,5,7,9,11.  The  remaining  36  receiver 
operating  characteristics  at  N  =  1024,  for  M  -  1,2,3,4,512,1024 
and  the  same  6  breakpoint  values,  are  presented  here.  Also,  a 
complete  set  of  characteristics  for  xQ  =  6  and  M  =  1,2,3,4,8,16, 
32,64,128,256,512,1024  has  been  added.  Finally,  a  couple  of 
isolated  examples  for  xQ  =  13  and  xQ  =  15  have  been  run. 

All  these  results  were  obtained  by  simulation;  the  number  of 
trials  utilized  is  indicated  on  each  figure.  A  couple  of  check 
cases  were  conducted  for  comparison  with  the  analytical  results 
in  [1].  They  are  xQ  =  3,  M  =  256  in  figure  B-ll,  which  verifies 
[1;  figure  39],  and  xQ  =  7,  M  =  8  in  figure  B-36,  which  verifies 
[1;  figure  11].  Labeling  on  all  plots  is  identical  to  figure  B-l. 

A  short  table  of  losses  relative  to  the  OPT-B,  for  N  =  1024, 
-3 

Pf  =  10  ,  P^  =  0.5,  is  presented  below.  These  were  extracted 

from  the  accompanying  receiver  operating  characteristics. 


xQ  dB 

1 

3 

5 

6 
7 
9 

11 

13 


loss  at  M  =  1 

5.5 

5.4 

4.2 

3.4 
2.7 

1.2 
0 
0 


dB  loss  at  M  =  1024 

0.5 

1.3 

3.1 

4.1 
5.0 
>7 
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Figure  B-3.  MGLR  ROC  for  xQ  =  1,  M  =  3,  N  =  1024 
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Figure  B-7.  MGLR  ROC  for  x  «  3,  M  *  1,  N  •  1024 


Figure  B-8.  MGLR  ROC  for  xo=3,M=2,  N=  1024 
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Figure  B-ll.  MGLR  ROC  for  x  =  3,  M  =  256,  N  =  1024 


67 


Figure  B-19.  MGLR  ROC  for  XQ  =  5,  M  «  1024,  N  =  1024 


Figure  B-23.  MGLR  ROC  for  xQ=6,  M=4,  N=  1024 
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Figure  B-27.  MGLR  ROC  for  xQ  =  6,  M  =  64,  N  =  1024 
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Figure  B-31.  MGLR  ROC  for  x  =  6,  M  =  1024,  N  =  1024 


Figure  B-35.  MGLR  ROC  for  x=7,  M=4,  N=  1024 
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Figure  B-47.  MGLR  ROC  for  x«ll,  M-4,  N>  1024 


Figure  B-48.  MGLR  ROC  for  x  =  11,  M  =  512,  N  =  1024 
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APPENDIX  C.  DERIVATION  OF  MODIFIED  LIKELIHOOD  RATIO  TEST 

For  known  (possibly  unequal)  signal  powers  {Sn),  the 
likelihood  ratio  may  be  shown  to  take  the  form 


LR  = 


1  +  *n 


1  +  — n 


Sn  2  0 


(C-l) 


leading  to  the  likelihood  ratio  test 


JL  sn  > 

S  1  +  S  xn  <  v 
n=  1  — n 


(C-2 ) 


However,  if  the  signal  powers  are  unknown,  we  replace  Sn  by 
hypothesized  value  Sn,  obtaining  the  modified  likelihood  ratio 

MtR  E  FT  |rr-r  rrV  xnll  -  sn  *  0  •  <c-3> 

n=l  (  n  n  ) 

This  quantity  is  maximized  by  the  choices  (random  variables) 
fx  -  1  for  x  ^  11 

Sn  *  \  ^  for  1  <  n  <  N  .  (C-4 ) 

(  0  for  xn  <  lj 

When  these  values  are  substituted  into  (C-3),  there  follows  the 
generalized  likelihood  ratio  test 

N  s 

YZZ  g(xn)  <  v  ,  (c-5 ) 

n=l  n  ' 

where  g(x)  is  given  by  (19).  This  is  a  familiar  result  [1]. 
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Under  hypothesis  Hq,  the  characteristic  function  of  z  is 
f°U)  -  fy(^)N  «  (Bj  +  —  B  *  exp(-b)  f  Bx  =  1-B .  (C-10 ) 

Upon  expansion,  there  follows 


f°K>  “  ZZ  (")  B f  »  Bn  - 1— - 

2  n=0  w  1  (1  -  i£)n 


(C-ll) 


The  corresponding  probability  density  function  of  z  under  HQ  is 


N 


n-1 


p°(u)  =  B^  S(u)  +  JZZ  (n)  Bi"n  fin  U-(-^i7?1~U)  for  u  *  °-  (C-12  ] 

Finally,  the  exceedance  distribution  function  of  z  under  Hq,  for 
threshold  v  >  0,  is 


Pr(z  >  v)  =  YZ.  (n)  Bi"n  Bn  Qn(v)  ' 
n=!  1  n 


( C— 1 3 ) 


where  functions  { Qn ( v ) }  are  defined  in  [7;  page  12, 
set  L  =  0  in  ( 23 )  ] . 

Expansion  (C-13)  is  a  very  useful  one  for  evaluation  of  the 
exceedance  distribution  function,  because  it  is  a  sum  of  positive 
terms,  of  which  the  components  are  quickly  and  easily  evaluated 
by  recursions.  Thus,  the  false  alarm  probability  can  be 
quickly  and  accurately  determined.  This  procedure  was  utilized 
for  the  determination  of  figure  25. 
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APPENDIX  D.  PROCESSOR  FORMS  FOR  VARIABLES  WITH  GAMMA  DENSITY 

In  this  appendix,  we  presume  that  all  the  bins  are  occupied 
when  signal  is  present;  that  is,  we  have  here,  M  =  N.  Therefore, 
independent  identically  distributed  random  variables  {xnJ, 

1  1  n  1  N,  have  common  probability  density  function 

a^  u^"1  exp(-au) 

P ( u )  =  - -  for  u  >  0  ,  (D-l) 

regardless  of  the  value  of  n,  where 

{a,  <  1  for  hypothesis  H..1 

f  ;  fJ  >  0  .  ( D-2 ) 

aQ  =  1  for  hypothesis  HqJ 

The  characteristic  function  and  cumulants  of  x  are 

n 

fxU) - - - -  ,  Xx(k)  =  (K-l)  1  for  k  >  1  ,  ( D-3 ) 

(1  -  H/a)^  x  ak 

respectively.  Thus,  Xx(l)  =  p/a*  and,  in  particular,  the  mean  of 
xn  is  p  (not  1)  for  hypothesis  Hq.  All  the  previous  results  in 
[1;  2;  3]  pertain  to  the  special  case  of  parameter  p  =  1  in 
(D-l).  (If  p  =  1/2,  I  integer,  then  random  variable  xfl  in  (D-l) 
and  (D-3)  is  a  scaled  chi-squared  variate  with  I  degrees  of 
freedom. ) 
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GENERALIZED  LIKELIHOOD  RATIO  PROCESSOR 


The  situation  is  identical  to  (D-l)  and  (D-2),  except  that 
strength- parameter  a^  is  unknown.  We  hypothesize  value  a^,  where 
a^  £  If  yielding  for  the  joint  probability  density  function  of 
observation  {xnJ  under  hypothesis 


—  P^ ( ' • 


/XN) 


N 

TT 

n=l  [ 


fal  xn_1  exP(“aixn) 


r(p) 


(  a 


P  nN 


r(p) 


n^-1  exp(-a1E), 


(D-6 ) 


where  random  variables 


N 

n  ■  77  x  , 

n=l  n 


N 

E  =  HZ  xn  • 

n=l 


( D-7 ) 


In  order  to  maximize  in  (D-6),  we  take  a  derivative  with 
respect  to  a-^,  obtaining 


tP-1 


r(P) 


N  af  1  exp(-aiE)  (pN  -  a^) 


( D-8 ) 


The  best  choice  of  a^  is  therefore  the  random  variable 


al 


if  E  >  pN 
1  if  E  <  pN, 


(D-9) 


because  a^  <  1  always.  The  resulting  maximum  of  p-^  in  (D-6)  is 
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because  useful  values  of  false  alarm  probability  P^  will  require 
threshold  values  v  in  (D-14)  larger  than  pN,  in  practice.  In 
fact,  the  mean  of  E  is  exactly  pN  under  hypothesis  Hq;  therefore, 
v  would  have  to  be  larger  than  the  mean  of  E  to  realize  small  P^. 

The  left-hand  side  of  the  processor  in  (D-14)  does  not  need 
or  utilize  ju  for  its  realization;  however,  achievement  of  a 
specified  P^  would  require  knowledge  of  fj  in  order  to  set  the 
correct  threshold  value  for  v. 

This  generalized  likelihood  ratio  test  is  identical  to  the 
optimum  processor  test  in  (D-5).  This  result  could  have  been 
anticipated,  because  the  likelihood  ratio  test  (D-5)  was 
independent  of  p  and  a^,  meaning  that  a^  was  never  needed  to 
realize  the  optimum  processor.  Again,  the  situation  here  is  that 
the  signal  occupies  all  bins  with  equal  powers  when  it  is 
present.  The  following  subsection  treats  the  case  of  unknown  and 
unequal  signal  powers  per  bin. 


95 


TR  10839 


(M/e)*' 
r Uj)  x 


for  n  e  L  , 


(D-18 ) 


n 


and  1  exp(-xn)/r(p)  otherwise.  Therefore, 


max 


px  =  TT 

neL 


r(p)  X 


n 


neL 


(p/e)p  exp(xn)  1 


n 


TT 

n^L 


N 

it 

n=l 


/xn_1  exP("xn) 


r  (p) 


xn-1  exP(~x«) 


n‘ 


T(p) 


( D-19 ) 


Meanwhile,  the  joint  probability  density  function  of  observation 
{xnJ  under  hypothesis  Hq  is,  from  (D-15)  and  ( D— 7 ) , 


_  n^-1  exp(-g) 

P°  r<„)N 


( D-20 ) 


The  modified  generalized  likelihood  ratio  follows  from  (D-19) 
and  (D-20)  as  (random  variable) 


MGLR  = 


max  p. 


(p/e)^  exp(xn) 


nsL 


n 


TTexp 

neL 


5s-  i 


-  In 


/■X  * 

H) 


=  expfp  YI1  g(xn/^)) 


(D-21 ) 


where  nonlinearity  g  is  independent  of  p  and  is  defined  by 


g(x)  = 


fx  -  1  -  ln(x)  for  x  >  bj 
[  0  for  x  <  bj 


( D-22 ) 
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